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The Einstein Telescope (ET) is a proposed third generation ground-based interferometric grav- 
itational wave detector, for which the target is a sensitivity that is a factor of ten better than 
Advanced LIGO and a frequency range that extends down to ~ 1Hz. Such a third generation in- 
terferometer will provide opportunities to test Einstein's theory of relativity in the strong field and 
will realize precision gravitational wave astronomy with a thousandfold increase in the expected 
number of events over the advanced ground-based detectors. A design study for ET is currently 
underway, so it is timely to assess the science that could be done with such an instrument. This 
paper is the first in a series that will carry out a detailed study of intermediate-mass-ratio inspirals 
(IMRIs) for ET. In the context of ET, an IMRI is the inspiral of a neutron star or stellar-mass black 
hole into an intermediate mass black hole (IMBH). In this paper we focus on the development of 
CNj ■ IMRI waveform models for circular and equatorial inspirals. We consider two approximations for 

the waveforms, which both incorporate the inspiral, merger and ringdown phases in a consistent 
^ | way. One approximation, valid for IMBHs of arbitrary spin, uses the transition model of Ori and 

rj^ Thorne lj to describe the merger and this is then matched smoothly onto a ringdown waveform. 

The second approximation uses the Effective One Body (EOB) approach to model the merger phase 
of the waveform and is valid for non-spinning IMBHs. In this paper, we use both waveform models to 
compute signal-to-noise ratios (SNRs) for IMRI sources detectable by ET. At a redshift of z = 1, we 
find typical SNRs for IMRI systems with masses 1.4M + 1OOM , 1OM + 1OOM , 1.4M©+500M Q 
q' 1 and 1OM +5OOM of ~ 10-25, ~ 40-80, ~ 3-15 and ~ 10-60 respectively. We also find that the 

two models make predictions for non-spinning inspirals that are consistent to about ten percent. 
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The Einstein Telescope (ET) is a proposed third generation ground-based gravitational wave (GW) detector for 
which the target is a sensitivity ten times better than that of the advanced detectors @, Q . By siting the instrument 
underground, it is also hoped that the range of frequencies to which the detector is sensitive can be extended into 
the 1Hz to 10Hz range, while also maintaining high frequency sensitivity up to 10 kHz. A design study for ET is 
currently underway, which is exploring the design, cost, site selection, etc., plus the potential scientific impact of 
such an instrument, with a view to maximizing the scientific output within a reasonable budget. If it is realized, 
ET will open up the possibility to study a wide variety of sources and address outstanding problems in fundamental 
physics, cosmology and astrophysics, e.g., studying compact objects (COs) and general relativistic instabilities, solving 
the enigma of gamma-ray-bursts and shedding light on their different classes, understanding the mass-spectrum of 
compact stars and their populations, and measuring cosmological parameters using GW sources as standard sirens 

The Advanced LIGO and Virgo detectors should achieve sensitivity in the low-frequency regime down to the seismic 
limit at roughly 10 Hz. For these detectors, the expected isotropic detection horizon (the angle-averaged distance 
to which a binary can be measured) for Advanced Virgo will be at ~ 150 Mpc (z = 0.035) for neutron star-neutron 
star (NS-NS) binaries and ~ 310 Mpc (z = 0.07) for neutron star-black hole (NS-BH) binaries ||. For Advanced 
LIGO, NS-NS signals will be in band for ~ 17 minutes and have a single detector detection horizon of ~ 200 Mpc 
(z = 0.045); NS-BH signals will be in band for ~ 4 minutes, and have a detection horizon of ~ 420 Mpc (z = 0.09). 
In contrast, if ET's seismic wall is pushed down to 5 Hz NS-NS signals will be in band for ~ 2 hours, and NS-BH 
signals for ~ 25 minutes. If the seismic wall is further reduced down to 3 Hz/1 Hz, then NS-NS signals will be in 
band for ~ 7hours/~ 5 days, and NS-BH signals will be in band for ~ 2 hours/~ 1 day, respectively. The Einstein 
Telescope detection horizon is expected to be at z ~ 1 for NS-NS signals, and at z ~ 2 for NS-BH signals 

The preceding paragraph served to illustrate that ET will be able to do the same type of science as the advanced 
detectors, but better. However, there is also some science which only ET will be able to do. The l-10Hz frequency 
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band lies below the range of current and advanced ground-based interferometers and above the range of the proposed 
space-based detector, LISA [6]. If ET achieves good sensitivity in this range, it will open up a new gravitational 
wave window. GWs in this range would be generated by sources with intermediate mass, i.e., from hundreds to a 
few thousand solar masses. Such intermediate-mass black holes (IMBHs) might be primordial, i.e., form in the early 
Universe as seeds from which BHs in galaxies subsequently grow, or they may form in the centre of dense globular 
clusters through runaway stellar collisions (see [3] and references therein for descriptions of these two scenarios). 
Observational evidence for the existence of IMBHs has accumulated over the last decade. This evidence is of two 
different types. First, ultraluminous X-ray sources (ULXs) have been observed that are not associated with active 
galactic nuclei and yet have fluxes many times the angle-averaged flux of a M < 20M Q BH accreting at the Eddington 
limit. Second, in several clusters, e.g., the globular clusters M15 and Gl, the stellar kinematics shows evidence for an 
excess of dark mass in the centre. In M15, there is also evidence for rotation in the core at a speed that is comparable 
to the central stellar velocity dispersion. One possible explanation for this rotation would be the transferal of angular 
momentum from a black hole binary to stars in the core through three body encounters. A binary comprising a 



20M Q black hole orbiting a 1O J M ob. 



ect at a semimajor axis of 10 pc would have enough angular momentum 
Ml for a discussion of this and other observational evidence for the existence 



to account for the observed rotation (see 
of IMBHs). 

The observational evidence will further improve with future X-ray and optical observations, but a robust mass 
determination of an IMBH candidate will be needed for a solid identification. This may be possible using radial 
velocity measurements of companions to ULXs in binaries, but this technique only yields a lower bound on the mass 
due to the unknown inclination of the system, and the companions are typically very faint. If IMBHs do exist, they 
will merge with other IMBHs and smaller COs, and these mergers will generate gravitational waves in the ET band. 
GW observations will provide very accurate mass measurements, and so the first convincing proof of existence for 
IMBHs may come from these observations Q. 

In this paper we will be concentrating on GWs from the mergers of IMBHs with lower mass COs, which might be 
either NSs or stellar- mass BHs. These sources are termed intermediate- mass-ratio inspirals (IMRIs). The dominant 
mechanism that leads to the formation of IMRIs is three-body hardening of an IMBH-CO binary in a core-collapsed 
globular cluster. 

Numerical simulations of globular clusters suggest that a fraction ~ 10~ 6 to ~ 10~ 4 of the ~ 10 6 initial stars that 
form a globular cluster will become stellar-mass BHs via normal stellar evolution Assuming a globular cluster 
with relaxation time of 1 G vr [lfj| . all these BHs should have formed within ~ 10 Myr, with the most massive BHs 
forming at around ~ 3 Myr 11]. These black holes should be more centrally concentrated than main sequence (MS) 
stars since there will be significant mass segregation of their higher mass progenitors [l2j], there will be preferential 
formation of stars near the cluster centre [13], and because BH birth kicks are not expected to displace BHs into 
the cluster halo [3]. Even if we assume that the BHs were distributed throughout the cluster, mass segregation 
should be able to assemble a subcluster of BHs near the centre after at most ~ 100 Myr. During mass segregation, 
BH-MS binaries will undergo three body and four body interactions that will replace the MS star by a heavier BH. 
Simulations [Tol. fl5jj suggest that, whether formed through successive BH mergers or stellar collisions, it is more likely 
to find an IMBH in a cluster with a dense core, so that IMBHs could exist in some tens of percent of current globular 
clusters. In this type of clusters, the first formation of an IMBH with mass ~ 100M Q could be^ 10 Myr after the 
subcluster of BHs is formed. IMBHs formed in this way have a negligible cross section for direct collisions/plunges 
with other objects. However, they will readily swap into a binary, as the heaviest object in the cluster, and such 
binaries can be driven to plunge through three body hardening or the secular Kozai resonance. IMBH-CO binaries 
will eventually be driven to merger through gravitational wave emission, and the GWs generated during these final 
phases of evolution could be detected by ground-based laser interferometers [7J. It has been estimated that in such 
environments the IMRI event rate for ET could be as high as a few hundred per year Q. 

We will describe two alternative models for the gravitational waveform generated during an IMRI, and then we 
will use them to estimate the signal-to-noise ratio (SNR) of events detected by ET. The low-frequency cut-off in the 
ET sensitivity will mean that only a short section of the inspiral phase will be observed for a typical ET IMRI event. 
A significant fraction of the total SNR will therefore come from the merger and ringdown and so it is important to 
include them in the waveform model. In both models, we will use the "numerical kludge" [l6| approach to describe 
the inspiral phase. This model was developed for extreme-mass-ratio systems for LISA and is based on the adiabatic 
inspiral of a test-mass through a sequence of geodesic orbits, with the sequence being determined by post-Newtonian 
energy and angular momentum fluxes, augmented with fits to Teukolsky fluxes and post-Newtonian conservative 



corrections [17] . The two models we consider here differ in the treatment of the merger. In the first model, which 



is valid for inspirals into IMBHs of arbitrary spin, we model the transition from adiabatic inspiral to plunge using 
the scheme provided by Ori and Thorne [l(. In the second model, which we use to cross-check our results in the 
non-spinning limit, we model the merger and plunge using the Effective One Body (EOB) approach. 

In our SNR calculations, we will assume that a network of ET-like detectors exists, rather than just a single detector. 
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This will be necessary if we are to estimate the extrinsic parameters of IMRI events to any precision. The duration 
of an IMRI event will be between a few seconds (e.g., a 10M Q plus 500M Q system) and a few minutes (e.g., a 1.4M 
plus IOOMq system). Over such a short timescale, the event is effectively a burst in the detector and so a single 
ET cannot determine the sky position of the source. While we do not consider parameter estimation in this paper, 
we assume the existence of a network since some of the science will require it. Paper II in this series will explore 
parameter estimation accuracies for the same types of ET network. 

This paper is organized as follows. In Section Ql] we describe the assumptions we have adopted for the detector and 
detector network, and describe the binary systems that will be used in the analysis. In Section [IIII we summarise the 
"numerical kludge" model that we shall use to model the inspiral phase. In Section HVl we summarize the transition 
to plunge scheme introduced by Ori and Thornc [1]. In Section [V] we describe the ringdown (RD) radiation and 
explain how the merger waveforms can be matched onto the RD. In Section PVTl we introduce the aspects of the EOB 
model relevant for our analysis. In Section IVIII we describe our implementation of the ET response function and 
describe the dynamics of twelve sample binary systems. We describe the orbital phase evolution up to the light ring 
(LR), i.e., the inner-most unstable circular orbit for massless particles, and present gravitational waveforms including 
the final inspiral, merger and RD. This section also includes a comparison between the two models for inspirals into 
non-spinning IMBHs. In Section IVIIII we will present ET SNRs for the sample binary systems and Section llXl contains 
our conclusions. 



II. ASSUMPTIONS 



A. Einstein Telescope Design 

In this sequence of papers, we take the response of a "single ET" to be that of two right-angle interferometers, 
coplanar and colocated, but rotated 45° with respect to each other. Assuming uncorrelated noise, this set-up is 
equivalent to the currently favoured ET design, i.e., a triangular configuration, with three 10 km long arms, and 
containing three independent detectors with 60° opening angles. This design has the capability to measure polarization 
at a single site, has lower infrastructure costs, and its sensitivity is a factor a/3/2v / 2 ps 1.06 higher than the two right- 
angle interferometer configuration [3] ■ We shall ignore this factor since it is small compared to other uncertainties in 
the design. 

In our studies, we will consider the "ET B" sensitivity curve, which at the current stage of the design study is the 
official sensitivity curve for ET [l8j . The corresponding amplitude spectrum is shown in Figure [T] 




FIG. 1: Sensitivity curve for the Einstein Telescope, as described in the text. The Advanced LIGO noise curve is also shown 
for reference. 



An analytic fit to the ET-B spectral density is given by [li 



aix bl + a 2 x b2 + a 3 x b3 + a 4 x b4 ] if / > f s 



(1) 



bh {f) -\oo if/</ s , 

where x — f/fa, f stands for the frequency, /o = 100Hz, Sq — 10 _50 Hz _1 , and f s is a low frequency cut-off that 
can be varied, and below which the sensitivity curve can be considered infinite for practical purposes. The various 
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coefficients take the values 

01 = 2.39 x 1CT 27 , h = -15.64, 

a 2 = 0.349 , 6 2 = -2.145, 

o 3 = 1.76 , 6 3 = -0.12, 

a 4 = 0.409, & 4 = 1.10. (2) 

Although it is hoped that ET will have sensitivity down to 1Hz, it is is not yet clear whether this will be achievable, 
and 3Hz might be more realistic. To be conservative, for the results described in this paper, we used a cut-off at 5Hz 
and started the inspiral evolution at the point when the GWs emitted by the IMRI swept through a frequency of 5Hz. 

As mentioned in the introduction, parameter estimation will require the existence of more than one, well-separated, 
detector. In this paper we will assume the existence of a detector network comprising three detectors sited at the 
current geographic locations of Virgo, Perth (Australia) and LIGO Livingston. We shall consider five configurations, 
C1-C5, for a few sample systems. These configurations are CI: one ET at the geographic location of Virgo; C2: as 
configuration CI plus a right-angle detector at the location of LIGO Livingston; C3: as configuration CI plus another 
ET at the location of LIGO Livingston; C4: as configuration C2 plus another right-angle detector in Perth; and C5: 
as configuration C3 plus another ET in Perth. 

In both this paper, and the next in the series, we shall first quote results for the most optimistic configuration C5, 
i.e., a network of three detectors each with the sensitivity of a single ET, and then we will go on to compare these with 
more modest configurations. The aim of this presentation is to exhibit the best performance achievable with a third 
generation ground-based detector network, but also to set the appropriate framework for the second paper, in which 
the assumption of a detector network will be required for the parameter estimation studies we will carry out. This 
presentation also follows the philosophy elsewhere in the literature, of first exploring the best performance achievable 
by detector networks, both in terms of detection and parameter estimation, for binary inspirals and burst-like events 
[19h21| . A 3 ET detector network is extremely optimistic, but we will show that more modest configurations will still 
produce fairly competitive results. 



B. Sample IMRI systems 



We will present results for twelve different binary systems. We take four combinations for the component masses, 
namely 1.4M© + 100M Q , 1.4M + 500M Q , 10M Q + 100M© and 10M Q + 500 Af , and three different values for the 
spin of the central IMBH, q — 0, 0.3, 0.9. The twelve sample binaries we examine are all possible combinations of 
these masses with these spin parameters. 



III. IMRI WAVEFORM MODELLING - INSPIRAL PHASE 



For comparable mass binary systems, the early inspiral phase is well modelled by post-Newtonian (PN) theory 
and the final few cycles can now be computed accurately using numerical relativity [22j . By contrast, extreme- 
mass-ratio inspiral (EMRI) systems, in which the mass ratio, r\ = 77111712/ (mi + TO2) 2 , is of the order of 10~ 5 , emit 
thousands of cycles in a regime where the velocity is a significant fraction of the speed of light. PN theory therefore 
does not apply, while numerical relativity cannot be used due to the large number of orbits that must be modelled. 
However, EMRIs can be accurately modelled using black hole perturbation theory, treating the mass ratio as a 
small expansion parameter. IMRIs lie somewhere between these two regimes, with mass ratios at which none of the 
preceding techniques have been tested. Accurate IMRI waveforms are therefore not known at present, and so we will 
construct our models using the best of what is currently available. We will take the inspiral model from the EMRI 
limit, but augment it with higher-order-in-mass-ratio PN corrections. We will consider two models for the merger 
and ringdown, and cross-check their predictions in order to improve our confidence in the results. While the resulting 
models may not represent the exact waveform of true IMRI sources, they will capture the main features of the signals 
and so should make reliable predictions for the SNRs and parameter accuracies that ET observations will achieve. 

We model the inspiral phase evolution, for both non-spinning and spinning IMBHs, using the "numerical kludge" 
waveform model described in [l7T |. This scheme was developed for the modelling of EMRI systems for LISA and 
has various nice features: a) the waveforms have been checked against more accurate, Teukolsky-based, waveforms 
for test-particles on geodesic orbits and the overlap exceeds 0.95 over a large portion of the parameter space 11611 : 
b) they are computationally inexpensive; c) conservative self-force corrections to this model have been derived [171 ] 
for Kerr circular equatorial orbits at 2PN order. This model is not complete as conservative corrections are not yet 
known for generic orbits, the phase space trajectories are approximate, although they have been matched to Teukolsky 
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based evolutions, and the waveform is constructed from the trajectory using a flat-spacetime wave-emission formula. 
Despite these various approximations, the numerical kludge waveforms should capture the main features of the inspiral 
waveform accurately. 

During the inspiral phase, radiation-reaction drives the motion of the CO. On short timescales, the small object 
follows an approximately geodesic orbit in the spacetime of the larger body. Over longer timescales, radiation- 
reaction causes the orbit to evolve adiabatically. This evolution can be characterized by changes in the geodesic 
orbital elements, namely the energy, E, angular momentum, L z , and the Carter constant, Q. In this analysis, we 
restrict our attention to circular-equatorial orbits, for which the energy and angular momentum take the form |23f 



E = 



E _ 1-2 (M/p) ± (o/M) (M/pf /2 
M v /l-3(M/p)±2(a/M) {M/pf 2 



L - z = _ ± / p y /2 1 T 2 (a/M ) (M/pf 2 + (a/M) 2 (M/pf ^ (3) 

mM \/l - 3 (M/p) ± 2 (a/M) (M/pf /2 

where p is the Boyer-Lindquist radius of the orbit, a is the spin parameter of the central black hole and the upper 
(lower) sign in Eq. ([3]) refers to prograde (retrograde) orbits. The Carter constant, Q, is zero for equatorial orbits. 

The orbital evolution of the CO is obtained by equating the rate of loss of energy E and angular momentum L z 
with the corresponding fluxes carried to infinity by the GWs, namely E and L z . These fluxes of energy and angular 
momentum must satisfy a consistency relation to ensure that circular orbits remain circular under radiation reaction 
fH, i.e., 



E(p) = ± / M L z (p) = n(p)L,(p), (4) 
p 6 i z ± oyM 

where d</>/di = H(p), is the azimuthal velocity of the orbit. The inspiral trajectory is given by 

p= dE E = dT z L - (5) 

since the energy and angular momentum fluxes are related to each other by (|4]), so we can use either the energy or 
the angular momentum flux to evolve the orbit. We choose to evolve L z , using a PN expression augmented by fits 
to more accurate fluxes c omp uted using the Teukolsky equation. The expression for L z takes the following form for 
circular equatorial orbits [25j 




336 \ p J \ p J 
high order Teukolsky fits | , (6) 

in which we have introduced q = a/M to denote the dimensionless black hole spin. The last ingredient to include in 
the inspiral phase evolution is the contribution of the conservative piece of the self-force (SF) to the evolution of the 
orbit. The conservative piece of the SF affects the frequency of an orbit at a given radius but does not lead to orbital 
evolution. The frequency shift accumulates over time, which affects the phasing of the waveform. We include this 
effect by modifying the <f> angular evolution as follows 

d^ _ /a> 

i.e., we multiply the phase derivative for a Kerr geodesic, labeled by the subscript "geo", by a frequency correction 
that depends on the instantaneous orbital parameters. The various conservative corrections contained in <5f2 can be 
derived by enforcing that asymptotic observables [26[ are consistent with post-Newtonian results in the weak field 
(27j . Once these corrections have been implemented the orbital frequency takes the form 17], 



1 + SQ), (7) 



dt p 3/2 ± ay fM \ ) 
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where the various coefficients are given by 

1 , 1975 , 27 , 191 , 1152343 , n . 

do= 8' dl= 896~' dl5 = ~10 7r ' ?1 - 5 = "160' d2 = 151584' (9) 

These are the ingredients that we use to build the inspiral part of our waveform model. In this paper, we will not 
explore the relative importance of the various terms that enter these expressions but leave that exercise for future 
work. 

Once the inspiral trajectory has been computed, the inspiral waveform can be obtained from an expansion of the 
form 

oo I 

h(t) = -(h+ - ih x ) hem -2Y em (8, $), (10) 

1=2 m=-£ 

where the spin-weight —2 spherical harmonics -2Ytm{0, are given in terms of the Wigner d functions by 



111 4- 1 

-,Y tm (0,*) = (-l) s y^dL(%™*, (11) 

k t (_i\k( ■ 0\2k+s-m( 6\2£+rn-s-2k 

with d l ms {6) = ^ {£ + m)\(£ - m)\{£ + s)\{£ - s)\ V jttj^ 2 j- 7n TTT, (12) 

vv ; ^ kl(£ + m - k)\(£ ~ s - k)\(s - m + k)l 



k=k, 



where ki — max(0, m — s) and kf = min(£ + m, £ — s). Additionally, the complex conjugates of the spin-weighted 
spherical harmonics satisfy 

a Y tm *(6,$) = (-l) s+m _ s y*- m (6»,$). (13) 

We include only the modes (£, m) = (2, ±2) in our model, which means that the components of the waveform at 
leading order are 

h+(t) = ^ [ ^ ' l f Jco B [2(0(t) + $)] J (14) 

K{t) = 4 M myp{t)] 2 cosgsin [2(0(t) + $)]; (15) 

where = dip/dt is the orbital frequency, p(t) is the radius of the orbit, and D is the distance to the source. 

This model provides a description of the inspiral only, but it ceases to be valid when the CO approaches the 
innermost stable orbit and adiabaticity begins to break down. We require a different model for the evolution of the 
CO through the inner-most stable circular orbit (ISCO) all the way down to the light ring, where we match on to 
ringdown radiation. 



IV. TRANSITION AND PLUNGE PHASES FOR AN INITIALLY SPINNING IMBH 

The approach we adopt to describe the merger waveform is based on the paper of Ori & Thorne [l| , which describes 
the transition from inspiral to plunge for circular-equatorial inspirals into a massive spinning black hole in the extreme- 
mass-ratio limit. To be consistent with their notation, we re-express the Boyer-Lindquist coordinates (t,p,9,(f>) in 
dimensionless form using the transformations r = p/M and t = t/M . 

As described in section IIII1 during the inspiral phase, the CO moves on a circular geodesic orbit with dimensionless 
angular velocity 

f2 = MO=§ = i 3 7 i— . (16) 
at r fi ' z + q 

As the CO inspirals onto the spinning IMBH, it radiates energy which is carried away by GWs. The radiation flux 
can be written as 

E GW = -E=^ V 2 n w / 3 £, (17) 
5 
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with £ being a general relativistic correction to the Newtonian, quadrupole-moment formula. 

The adiabatic inspiral phase evolution continues until the CO approaches the ISCO. As we discussed in the previous 
section, the adiabatic prescription breaks down somewhat before ISCO, as the orbit starts to evolve more quickly 
and the instantaneous-geodesic approximation is no longer valid. Hence, we need to find some point near the ISCO, 
Straus <; HSCOj that joins smoothly the adiabatic inspiral of Section IlIII onto the transition phase. This choice must 
ensure a continuous and smooth matching of the inspiral and transition waveforms. The transition radius, ftrans ; 
is the point at which ~ri nsp irai becomes "too fast" . "Too fast" obviously has a different meaning according to the 
binary system under consideration as it depends on the mass-ratio r\. However, the transition solution described in 
the following provides a unique solution that should match smoothly onto the inspiral. In practice, we chose a fixed 
matching radius for each binary system in order to give a smooth matching onto the transition waveform, but have 
subsequently verified that changing the value of~ri nsp ; ra i at which the matching is done does not significantly alter the 
resulting waveform, i.e., our choice of ft rans is robust. 

At the transition point, the circular geodesic has dimensionless angular velocity, energy and angular momentum 
given by 

n tIans = MCl= , (18) 

trans ' Q 



£S -Etrans -Etrans 1 2/ft rans + ^/^trans ;, n \ 

-^trans = — 7T — — , , (T y J 



\J 1 - 3/ftrans + 2q/r^ ns 



^trans _ -1/2 1 ~ V^t/aL + ^ l f trans / 9n N 

M Vl-3/f teanB + 2 g /f^ 2 



rans 



The values of these quantities for the binary systems under consideration are given in Table HI We have included the 
values of these quantities evaluated at ISCO for reference. 



Binary systems 


q = 0.9 


g = 0.3 


L z 


E 


'trans 


U 


E 


'trans 


[10+100] M Q 


2.106 


0.8477 


2.450 


3.157 


0.9309 


5.124 


[1.4+100] M Q 


2.104 


0.8472 


2.442 


3.155 


0.9308 


5.120 


ISCO 


2.100 


0.8442 


2.321 


3.154 


0.9306 


4.979 


[10+500] M Q 


2.104 


0.8472 


2.428 


3.155 


0.9308 


5.075 


[1.4+500] M Q 


2.104 


0.8472 


2.427 


3.155 


0.9308 


5.075 



TABLE I: Dimensionless values for the energy E and angular momentum L z as denned in Eq.Q at the point of transition 
ftrans. The values for the energy and angular momentum at ISCO have been included for reference. 



As the CO enters the transition regime, the geodesic motion ceases to be adiabatic but radiation reaction continues 
to drive the orbital evolution. In this regime the CO still moves on a nearly circular orbit with radius close to ftrans- 
Additionally, since the radiation reaction is proportional to the mass ratio and is therefore weak, the angular velocity 
and proper time can be approximated by |l| 

f2 ~ ftrans , (21) 



dt 



-. 



df / df \ V 1 - 3/ftrans + 2q/f^ ns 



(22) 



dt ~ \dtj trans l + q/?HL 

In the vicinity of f tra ns) the CO's energy and angular momentum can be written 

E = Et rans + fitransf j L = -i/trans + £, (23) 



8 



where 

^ = — kt\ , and (24) 
dr 

k - —n 7/3 1 + i/^L 6 (25) 

§ "trans i ''trans ■ l/OJ 

Y 1 - 3/ f trans + 2?/ r t / ans 

To make further progress, we recast the effective potential describing the radial motion for geodesies 

V(f,E,L) = E 2 ~^([E(f 2 + q 2 )-L q ] 2 -(f 2 -2f + q 2 )[f 2 + (L-Eq) 2 }) , (26) 

as a function of f and £ = L — Ltrans- We use the variable R = f — ft rans to parametrise the CO's location during the 
transition regime. Both R and £ are small and therefore we can expand the effective potential in terms of these two 
variables 

V(R, = ^R 3 - 2PR£ + constant , (27) 

where a and j3 are constants to be computed below. This new expression for the effective potential can be plugged 
into the equation describing the radial evolution in the transition regime, namely, 

d 2 r l dV(r,£) - 

d^--2—d^ +VFseW - (28) 

The radial self-force r)F sc \f is approximately non-dissipative and hence can be ignored. Absorbing this term into 
— \dVjdr effectively amounts to changes in r trans , -Et xans , -^trans and a by fractional amounts proportional to r\. 
Likewise, these various quantities change by order 0(77) as a result of the CO's perturbation of the black hole's 
spacetime geometry 28]. We shall ignore these small corrections in the following analysis. 

Once we insert the effective potential ([27| into the equation of motion (|28|) and set f = at the moment when 
£ = 0, we obtain [l[ 

^ = - a R 2 - V (3ht , (29) 



where the constants a , j3 are given by 



1 ( d 3 V(r,E,L) \ 

a= l[ ' (30) 



_1 *Vi?,E,L) (31) 
2 I dLdf dEdf ) 

\ / trans 

Eventually the transition regime ends, radiation reaction is no longer important and pure plunge takes over. There- 
after, the CO plunges towards the black hole with nearly constant energy and angular momentum given by [l[ 

ifin — ^trans = — (^ToTplungc)?? 4 ^ 5 , 

Sfin — £trans = — ^trans(KToTpi ungc )?7 4 / 5 , (32) 

where, 

Tpi ullg0 = 3.412 , T = (a/3 K )~ 1/5 • (33) 
During the plunge phase, the evolution is given by the Kerr geodesic equations [23], i.e., 

fr_ _ 6 g fin L fin q + L 2 n (f - 3) + (g 2 - rjr - E 2 n g 2 (r + 3) 

df 2 " f 4 > W 

d^ L na (f-2) + 2E Rn q 

(35) 



dt ' E Rn {P + (2 + f)q 2 )-2qL Rn : 
df r(q 2 +r(r- 2)) 

di ~ E Rn (P + (2 + r) q i)-2qL Rn ' 



(36) 
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Notice that the plunge angular frequency (|35[) is entirely determined by the values of the energy and angular momentum 
Q32[) . Hence, to match the transition regime onto the plunge regime, we only need to find the point fpi unge at which 
the transition angular frequency (fTB"|) and the plunge angular frequency (|35|) smoothly match for these specific values 
of energy and angular momentum (|32l) . 

We now have all the ingredients required to build our waveform model from inspiral to plunge. For the inspiral 
waveform the cross and plus polarizations are given by (p"4"l) . (fl"5j) . which provide a good approximation before ftrans- 
Thereafter the particle is no longer on a circular orbit. Hence, we shall use the following expressions 



1 1 - 2 cos 26 cos 2 [<f>(t) + $] - 3 cos[2(<£(i) + 



(3 + cos26»)|2cos[2(0(i) + $)]</> 2 '(i) + sin[2(</»(t) + $)]</>(/;) }r 2 + |4(3 + cos 2(9) Bin [2 (<£(*) + 
(1 - 2 cos 26* cos 2 + $] - 3 cos[2(^(i) + $)])rj? 

— 2/Lt COS 6* 



(37) 



D 



sin[2(0(i) + $)]r 2 + { cos[2(0(t) + - 2 sin[2(</>(i) + $)]</> 2 '(i)}? 



+ |4cos[2(0(i) + $)]<£(i)r + sin [2 (<£(*) + $)]f }r 



(38) 



where zx is the mass of the CO and D is the distance to the source. These are flat-spacetime emission formulae applied 
to a geodesic in curved space and so is in keeping with the "numerical kludge" approach to waveform generation. This 
approximation is vindicated by the similarity of the results to the EOB waveforms which we will demonstrate later. 

This "transition" waveform model provides a consistent modelling for the gravitational radiation emitted from the 
early stages of inspiral evolution all the way down to the horizon. However, we shall attach the final part of the 
waveform, i.e., the RD part, at the effective light ring. The following section describes the method for attaching the 
plunge phase onto a RD waveform. The method is generic, i.e, it is applicable both for "transition" model, and for 
the EOB scheme which we will describe in Section IVII 



V. RINGDOWN WAVEFORM 



The ringdown radiation originates from the distorted Kerr black hole that is the end product of the merger, and 
consists of a superposition of quasinormal modes (QNMs), labelled by indices (£,m,n), where (£,m) specifies the 
mode and n the tone. Each mode has a complex frequency ui, whose real part is the oscillation frequency and whose 
imaginary part is the inverse of the damping time [30J , 

L) = LOlmn — i/Tlmn- (39) 

These two observables arc uniquely determined by the mass and angular momentum of the newly formed Kerr black 
hole. Recent numerical simulations have shown that the total mass-energy radiated during the merger of two equal- 
mass maximally spinning BHs ranges from 0.6% — 5% of the total rest mass energy. The energy released in IMRI 
ringdown radiation will be much lower, as it is suppressed by the mass ratio r\. 

It would be reasonable to ignore the change in mass of the IMBH and approximate the final mass and spin by the 
initial values for the central black hole. But, for completeness, we will use a one-parameter fit, derived within the 
framework of the EOB formalism, as an approximation for the final mass Mt of the black hole even though this does 
not account for the IMBH spin [311 ] . 



M f /M = 1 + (^/8/9 - l)r] - 0.498(±0.027) if. (40) 

This fit is consistent with NR simulations to about ~ 2% accuracy for mass ratios r\ > 0.16 [Hj]. However, the 
extrapolation to smaller values of 77 is consistent with NR simulations [H| and test-mass limit predictions [32] • Note 
that the coefficient of the linear term has been fixed to the test-mass limit value. 

We will also use the fit by Rezzolla, et. al. [13], to compute the final spin of the black hole 

a f /M f = qj = q + s 4 q 2 r) + s 5 qr) 2 + t [) qi 1 + 2 VSr] + t 2 if + t 3 if, (41) 

in which the coefficients, obtained via a least-squares fit to available data, are 



s 4 = -0.129 ± 0.012, s 5 =0.384 ±0.261, 

t = -2.686 ± 0.065, t 2 = -3.454 ± 0.132, t 3 = 2.353 ± 0.548. 



(42) 
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This fit was derived by using available data for spin parameters q <; 0.8 and mass ratios 77 ^ 0.16 along with exact 
results which hold in the extreme-mass-ratio limit, i.e., 77 — >■ 0. For q — 1 the fit (|41[) is a non-monotonic function 
with maximum qt — 1.029 for rj ~ 0.093, but this fit should be valid for the systems we consider here, for which 
i] <^ 0.08 and q < 0.9. At present we do not have a good idea as to plausible values of q for the systems that will be 
detected by ET. If the dominant process by which IMBHs acquire spin is through the capture of COs, their angular 
momenta will undergo a damped random walk (35l |. [36j]. This process was studied in detail by Mandel [37'], who 
computed the probability distribution for the spin of IMBHs that gain mass following a series of minor mergers. This 
work suggested that IMBHs with masses in the ~ 10 2 -10 4 M Q range would have spin parameter q ~ 0.3. 

Given the mass and spin of the final Kerr BH, we can uniquely determine the complex ringdown frequencies (|39|) . 
We do this by building an interpolation function based on the data provided in Table 2 of Berti, et. al. [3a |. and 
evaluating this function for the final IMBH spin. Following Berti, et. al. [Hj], and Buonanno, et. al. |30|, we construct 
a RD waveform that includes the fundamental mode (£ — 2, to = 2, n = 0) and two overtones (n = 1, 2). In order to be 
consistent, we also include the "twin" modes with frequency w' imn — —W£- mn and a different damping T' tmn — r^_ mn . 
We include the "twin modes" since a mode with a given (£, to) will always consist of a superposition of two different 
exponentials. It might be the case that one of the exponentials has a shorter damping time or is less excited in 
the given physical situation and hence become "invisible" , but formally we cannot have an isolated u £ = m = 2" 
frequency with a positive real part. Furthermore, a single-mode expansion restricts attention to circularly polarized 
GWs. These considerations are not important for non-spinning BHs as the two mirror solutions are then degenerate 
in the modulus of the frequency and in the damping time. Adding black hole rotation acts in a similar way to an 
external magnetic field on the energy levels of an atom, causing a Zeeman splitting effect of the QNM frequencies. 
The ringdown waveform is given by 



h(t) = -{h+ - i 



Iran 

+ ^ m ne i( """" t+ ^'" J e- t/Tto '"5; m (a^7 mn ) 



(43) 



in which Mf and D are the mass of the Kerr BH formed after merger and the distance to the source, respectively. 
The spherical harmonics of spin-weight —2, - 2 Si m obey the equation 



1 d 



9 9 . (to — 2 cos 9Y „ 
aw sin 9 - — + Aaw cos 8 - 2 + 2maw + A 



. sin 9 d6 1. av J sin 
Expanding _2S^ m (aw) and the eigenvalue A [39I ] 

-iSi m {aw) = ^Yem + awS^l + iaw^S^ + O^aw) 3 ), 
A = A + aw\ 1 + a 2 w 2 \ 2 + 0((awf) 7 



- 2 S tm (aw) = 0. (44) 



(45) 



where -%Yt m are the spherical harmonics of spin weight s = —2. The normalizations of -lYim and - 2 Si m {aiw) are 
fixed by 



I - 2 Ye m I ^sin 9d9 = / \- 2 S em (aw)\' sin 6 d6 = 1. 
^0 



(46) 



After plugging (|45j) into (j44|) and collecting terms, we obtain the zero and first order corrections to the eigenvalue A, 
i.e., 



A = (£- l)(i + 2), Ai = -2m 
The first order correction to - 2 Sg m (aw) is given by 

_ \^ J! Vl 



1(1+1) +4 

£[£ + !) 



I'm , 



(47) 



where the non-zero coefficients <r„ m are 



(e + iy 



(£ + 3)(l - l)(t + m+ !)(£ - to + 1) 



(2£+ l)(2£ + 3) 



1/2 
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(l + 2)(l-2)(l + m)(l-m) ni/2 

(2£+l)(2£-l) J ' 1 ' 

This first order approximation will suffice for our analysis since auj < 1 in Eq. (1441) for the systems we shall 
consider. Using these relations we can split (|43p into plus and cross polarizations, which can be matched onto the 
corresponding components of the plunge waveforms. This matching requires the determination of 24 constants, 12 for 
each polarization, but the equations can be rewritten as two 6D systems of equations. We do this matching at the 
time that the particle reaches the light-ring. 

The approach to determine the various amplitudes and RD phases of Eq. (|4"5|) is the following: use the plunge 
waveform to compute ten points before and after the LR to build an interpolation function. This function can be 
used to match onto the various QNMs by imposing continuity of the waveform and all the necessary higher order 
time derivatives. We first match the plunge waveform onto the leading RD tone n — at the point where the orbital 
frequency f (|3"5j) . ([55]) ) peaks, i pC ak- This fixes 2 constants per polarization. We then use these values as a seed 
to compute the amplitudes and phases of the first overtone at a time t poa k + dt. Finally, we use the values of the 
amplitudes and phases of the leading tone and first overtone to determine the four remaining constants at a time 

^pcak ~\~ 2tt£. 

In terms of the relative SNRs contributed by the different tones, we have found that the leading tone along with its 
twin mode are the major contributors. These two modes contribute more than 90% of the SNR for the various binary 
systems under consideration. The remaining two overtones and their twin modes do not contribute substantially to 
SNR but we have included them for completeness and to achieve the best possible modelling of the RD waveform. 



VI. PLUNGE AND MERGER WAVEFORM FROM THE EFFECTIVE ONE-BODY APPROACH 

In order to cross-check the predictions of our IMRI "transition model" in the non-spinning limit, we have used 
an independent approach to model the plunge and merger phases, which is based on the Effective-One-Body (EOB) 
formalism. Note, however, that the the ringdown waveform presented in Section [V] is generic and can be applied to 
this new prescription. Therefore, within this framework, a complete waveform model consists of the inspiral waveform 
described in Section lHTl the plunge and merger waveform to be discussed below, and the ringdown waveform described 
in Section IVl 

Although we use the EOB model only in the non-spinning limit, an extension of the EOB scheme does exist 
which includes leading-order spin-orbit and spin-spin dynamical effects of a binary system for an "effective test 
particle" moving in a Kerr-type metric [ibj . and next-to-leading-order spin-orbit couplings [4l[. However, it has 
been recently found [52] that it is not straightforward to include higher-order non-spinning PN couplings, such as the 
4PN and 5PN adjustable parameters that were recently calibrated to numerical relativity simulations for non-spinning 
systems [HI, [Hj] , using these Hamiltonians [io|, |4l[ . Additionally, the EOB Hamiltonian in fjljdoes not reduce to the 
Hamiltonian of a spinning test particle in Kerr spacetime. This issue was recently resolved in [45| . in which a canonical 
Hamiltonian was derived for a spinning test particle in a generic curved spacetime at linear order in the particle spin. 
The construction of an improved EOB Hamiltonian based on the results of [45[ has recently been obtained 46] . The 
Hamiltonian derived in [41| has recently been used in an exploratory study to calibrate the EOB parameters using 
numerical relativity simulations of spinning, non-precessing, equal mass BHs. This is the same approach that was 
previously used with great success for non-spinning black hole systems, e.g., to derive fits for the final mass and spin 
of a BH after merger that are consistent with NR to about ~ 2% accuracy [Hj]. We used the non-spinning EOB 
model only in this work because that model is more mature. However, the EOB model has recently been used to 
model circular-equatorial extreme-mass-ratio inspirals (EMRIs) around spinning supermassive black holes, using fits 
of various post-Newtonian parameters to Teukolsky-based waveforms [471] ■ Comparisons between an IMRI model 
based on this spinning EOB framework and the waveforms constructed using the transition model should be pursued 
in the future. 

The standard analytic method to study the two-body dynamics of comparable-mass BHs is the PN framework, 
which is an expansion in the characteristic orbital velocity v/c. Accurate equations of motion at 2.5PN order were 
derived in the 1980's (48[. At present, corrections at 3.5PN-level in the equations of motion are available [49l - [5lj ] and 
the diffeomorphism-invariant dimensional regularization method proposed by Blanchet, et al., [52J opened up the way 
to also study the GWs emitted by inspiralling non-spinning compact binaries up to 3.5PN order of accuracy. The 
PN formalism is reliable as long as the PN expansion parameter d = M/d -C 1/6 , where M is the total mass of the 
binary system and d is the separation between the two BHs, i.e., only during the early inspiral stage. Once S ;> 1/12, 
we require an alternative description of the motion and radiation to accurately model the final stages of inspiral plus 
merger and ringdown [22j ]. Numerical Relativity (NR) is the best candidate to model these various evolution phases 
for comparable-mass BHs. However, the need to compute thousands of waveform templates to carry out matched 
filtering searches makes NR an impractical tool due to the high computational cost of producing individual waveforms. 
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In order to circumvent this problem, Buonanno & Damour |53| introduced a framework to study the entire waveform 
using conservative dynamics at 2PN order. This scheme provides an accurate model from the early stages of inspiral 
down to the LSO. Inside the LSO, the scheme includes a non-perturbative plunge down to the LR and, thereafter, 
the model can be matched onto ringdown radiation. The basic claim is that it is possible to use analytical tools to 
obtain a sufficiently accurate waveform from inspiral to RD. Instead of using PN expansions in the original form, 
resummation methods (a non-polynomial function of the mass ratio r/ = rn\m2/M 2 , which incorporates some of the 
expected non-perturbative features of the exact results) can be used to improve the convergence properties of the PN 
expansions. In the test-mass limit (77 — > 0) the two basic ingredients of a GW signal, the two-body energy and the 
GW energy flux, are first resummed. The resulting EOB model is constructed to exactly recover geodesic motion in 
the test-mass limit. The argument was that the resummed quantities would also provide a good description for the 
comparable-mass case, since this is in effect a smooth deformation of the test-mass limit. The intrinsic flexibility of 
the model allowed a natural extension to higher order [54| once the 3PN calculation was completed. This belief has 
been vindicated by recent advances in numerical relativity, as it has proven possible to construct EOB waveforms 
that match the results of the numerical simulations very well [22j. For instance, the EOB model is the only analytic 
approach that has successfully predicted the spin of the final BH after a merger to better than ~ 2% accuracy |32j . 

In building the model, one of the challenges was to encode the conservative part of the relative orbital dynamics 
into the dynamics of an effective particle, i.e., to map the real conservative two-body dynamics at the highest PN 
order available, onto an effective one-body problem, in which a test particle of mass /1 = mxmij M moves in some 
effective background metric gf^,. In the particular case of non-spinning binaries and ignoring radiation-reaction 
effects, the best effective metric was found to be a deformation of the Schwarzschild metric, with 77 playing the role 
of a deformation parameter. 

We shall now introduce the mathematical machinery required to construct the plunge waveform model. We first 
show how to obtain the cross and plus polarizations for the plunge phase evolution and then explain how we match 
the inspiral waveform model of Section HTT1 onto this plunge waveform. We shall use phase space variables (r, <fi,p r ,P<i>) 
to write the effective EOB Hamiltonian for non-spinning binaries at 3PN order as follows 



H eff (r, p) = fj,H eS (r,p) =fj,* A(r) 



1 + p 2 + 



( Mr) 
\D(r) 



1 (n-p) 



(2:1 (p 2 ) 2 + z 2 p 2 (n • p) 2 + z 3 (n • p) 4 ) 



(49) 

where n = r/r and r = |r|. It is convenient to replace the radial momentum p r by the momentum conjugate to 
the tortoise radial coordinate 7% = J dr (_B(r)/A(r)) 1 ^ 2 , where A{r), B(r) are metric functions which will be defined 
below. We do this because p rt tends to a finite value after merger, whereas p r diverges at the event horizon. Such a 
coordinate transformation allows a more controlled treatment of the late part of the EOB dynamics [22| • Under this 
coordinate transformation the EOB Hamiltonian takes the form 



\ 



Pr, 



and the mapping between the real and effective EOB Hamiltonian is given by 

1 



Hvon{r,Pr,,4>) = 



rrrcal 
-"EOB 



n 



l + 2r7(iJ cff -l), 



(50) 



(51) 



where -ff roa i = H TC!L i/fi. This relation holds true at 2PN and 3PN order. The arbitrary coefficients 2:1,2:2 and 2:3 in 
Eq. (|49|) are subject to the constraint 



821 + 4z 2 + 3z 3 = 6(4 — 377) 77 . 



(52) 



It is possible to determine these coefficients by means of a fit to numerical results for comparable mass systems. We 
should bear in mind that in so doing, we must recover exact results in the test-mass limit. To achieve this, we have 
to ensure that z%, 22 and 23 must go to zero as rj — > 0. On the other hand, Damour, et al., [55j |. found that the terms 
proportional to 22, 23 in Eq. (|49[) are very small for quasicircular orbits. They also noticed the convenience of setting 
z\ = because in any other case Z\ could be suitably chosen so as to cancel the 3PN contribution in the metric 
coefficient A(r). Hence, we shall follow the general philosophy adopted by Damour et. al. |55j . and Buonanno et. al., 
[5(| and set Z\ = z 2 = 0, z 3 = 2(4 — 3?/)?/. In order to ensure the existence and 77-continuity of a last stable orbit as 
well as the existence and 77-continuity of an ^-deformed analog of the light-ring, the metric coefficient A{r) must be 
Pade resummed, 
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Af^(r) = 



where 



r 2 [(04(77) + 877 - 16) + r (8 - 277)] 



r 3 (8 - 277) + r 2 [04(77) + 477] + r [204(77) + 877] + 4[?7 2 + 04(77)] 



(53) 



04(77) 



94 41 , 
¥-32^ 



(54) 



The light-ring in the test-mass limit is the solution to d/dr(A(r)/r 2 — 0). For 77 ^ 0, the 77-deformed LR is obtained 
by solving d/dr(A(r,r])/r 2 ) = 0. The existence of this "deformed" LR guarantees that in its vicinity the orbital 
frequency ft reaches a maximum. Additionally, the Pade resummation of the metric coefficient D(r) ensures an 
77-continuity in the plunging phase. Its Pade resummed form at 3PN is given by 



6rjr + 2t?(26 - 3t?) 



(55) 



Using reduced quantities H Ica i = H tea \/[x, t = 
waveform evolution from the LSO to the LR are 



Q.M, the EOB equations of motion that describe the 



d0 
"dl 

dr 
d7 

di 

di 



Apt, 
i]r 2 HH cS 

A^ 1/2 
B 



1 



r)HH cS 



2A , 



— J* A 



1/2 



1 



2 V HH cM 



A' 



A' 



t) 



Z3 



A' 2A\ 



(56) 

(57) 
(58) 
(59) 



where J-^ is the <fi component of the radiation-reaction force. 

We shall use this set of equations only to model the waveform evolution from the LSO to the LR and it is therefore 
reasonable to assume that the two-body dynamics is no longer driven by radiation-reaction but occurs along a geodesic 
with constant angular momentum given by 31 j 



A'(u) 



(u 2 A(u))' 



(60) 



LSO 



where u = 1/r and the prime denotes d/du. The motion during plunge remains quasi-circular in the sense that p 2 
stays numerically small compared to pi. 

In the limit rj 0, we know that circular orbits in a Schwarzschild geometry satisfy the Kepler law Q 2 r 3 = 1. It is 
common to use = Q 1 / 3 or xq = Q 2 / 3 to describe all PN corrections whether they are proportional to the square 
of the linear azimuthal velocity (v$ = fir) or to the gravitational potential (U = 1/r). In order to generalise this to 
the case rj ^ 0, we introduce the functions 



dA{r) 
dr 

ril){r,p<t,), 



1 + 2,/ I , Mr) { I '-}-!] • (61) 

'•;> : : r c( r.p,,). (62) 

so that this modified radius t-q is related to fi by the standard Kepler law fi 2 r 3 1 = 1. This expression holds true during 
plunge, while the combination K — tt 2 r 3 becomes of order 0.5 at the effective LR. We use this modified relation, 
1, to determine fi inside the LSO. Finally, we use Eqs. (4.38) and (4.39) of [57j to determine the value of p Tt 



2^3 _ 



at the LSO from which we can solve the set of Eqs. (|56|) . (|57|) and ([59)) for the plunge phase. 
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Once the orbital evolution has been determined, the EOB waveform is given in terms of spin- weight —2 spherical 
harmonics _2^ m (^, 3>) through the relation 

htm = ~(h+ - ih x ) em = - J dft - 2 Y e * m (8, $) (h + - ih H ) . (63) 

As for the inspiral, we shall include only the two modes (£, m) = (2, ±2). At leading PN order the two modes h 22 and 
/12-2 are given by 

^22° B « = -S^vim^m) 2 F 22 (t)e- 2i +M , (64) 

h^° 2 B (t) = h* 22 (t) (65) 

where D is the distance to the source, M is the total mass of the binary system and 4>(t) is the binary orbital phase. 
The factor F 22 is a resummed version of all the PN corrections and is given by 

F22 (t) = H cS T 22 {t) p 2 22 {x{t)) e lS ^ l \ (66) 

where x(t) = r^ 1 and Tg m (t) is a resummed version of an infinite number of logarithmic terms that enter the transfer 
function between the near-zone and far-zone waveforms. These terms arise due to tail effects connected to the wave 
propagation in a Schwarzschild background of mass -Madm = Heob 1^3 ■ The factor Si m is a supplementary phase 
which corrects the phase effects not included in the complex tail factor Ti m . 

Finally, in order to enhance the agreement between the EOB model and numerically computed waveforms near the 
end of inspiral and during the beginning of plunge, we have introduced the resummed quantity p£ m , which enters 
the waveform only through its I- th power, p\ m . Previous waveform models utilised a different PN improving factor 

F 22 (see Eq. ([66]) ). namely, F 22 (t) = H cS T 22 (t) f 22 (x(t)) e iS22 ^\ where / 22 is a PN-expanded amplitude factor. The 
Taylor-expanded fim's produce results that are incompatible with numerical data close to the LSO. This problem 
arises because the /^ m 's have coefficients that grow linearly with £, and these terms are problematic for the accuracy 
of PN-expansions, as shown in 58]. Replacing ft m by its g-th root pi m — [ftm] 1 ^ seems to be a cure for these 
accuracy problems and improves agreement with Numerical Relativity in the strong-field/fast-motion regime [22j . 
The explicit forms of the various quantities introduced above are [58l . l59j , 

T 22 = r(3 r 7 3 f fc) e7rle2illos(2fcro) ; (67) 

622 = -H Ieal Q+—Tr(H leai nf-2Ar 1 x 5 / 2 , (68) 
d 105 



. . /55f 43 \ /19583i/ 2 33025i/ 20555\ 

= i+^—--j x +^-^—-——--— l j x 



2 



/10620745^ 3 6292061^ 2 41tt 2 ^ 48993925^ 428 1556919113\ 3 

V 39118464 3259872 + ~~L92 9779616 ~ 105 CU ° r ° S2 ^ + 122245200 J X 

/9202 , , . . 387216563023\ 4 / 439877 , , , . 16094530514677\ 5 m , 6 , 

culerlog,(a;) \x + eulerlog, (x) x° + 0(x b ), 

\2205 &2V ' 160190110080/ \ 55566 533967033600 / v ; ' 



(69) 



in which 



k — GH rea \mQ, k = mQ, 
r Q = 2M, x(t) = — , 



1 



eulerlog 2 (a;) = -y E + 2 log 2 + - logx, with j E = 0.577215, 



(70) 



with M being the total mass of the binary as before. 

In order to facilitate the matching between the inspiral and plunge waveforms we rewrite the factor F 22 (t) 
G{t) e le w, where both G(t) and e(t) are real quantities. The waveform can then be written as 
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h(t) — —(h+ — ihx) 

= - 2 Y 22 (6, $) h 22 (t) + - 2 Y 2 _ 2 (8, $) h 2 _ 2 (t), 

and Eq. ()T1[) allows us to write the cross and plus waveform components for plunge and merger as 



h+(t) 



AjjLG{t)[Q.{t)r n (t)f /1 + cos 2 



D 



2(0(f)+*--e(f)) 



4/iG(t) r n (f)] 2 . 
/i x (C = — cos 6* sin 



2 + * ~ 2 eW) 



(71) 

(72) 
(73) 



Now that we have found explicit expressions for the plus and cross polarizations of the plunge waveform, we must 
match these on to the inspiral waveform, which is achieved as follows: we already know that the scheme used to 
evolve the radial coordinate in the inspiral phase, Eq. ([5]), begins to break down as the CO nears the ISCO. Hence, 
we first need study the behaviour of Eq. ([5]) near the ISCO to find out where the adiabatic approximation starts to 
break down. We have found that for the binary systems under consideration, the transition to plunge occurs at a 
point rt rans ~ Hsco ■ In practice, the matching point /Vans was chosen to ensure that: i) the transition from inspiral 
to plunge of the radial, Eqs. ([5]), (f57| . and azimuthal, Eqs. (j8|), (f56|) . coordinates is smooth; ii) the transition from 
inspiral, Eqs. (fl4"[) . (|72|) . to plunge, Eqs. (|15l) . (f73"[) . of the waveform is smooth; and iii) the choice of r trans is robust, 
i.e., conditions i) and ii) are met in a vicinity around the precise value of rtrans- Furthermore, we have verified that the 
phasing and the amplitude of the resulting waveform are insensitive to the exact choice of rt ran s- Table HTl indicates 
the energy and angular momentum, as defined in Eq. ([3]), at the transition point r trans compared to the corresponding 
values for a test-particle at ISCO. 



Binary systems 


Lz 


E 


^"trans 


[10+100] M @ 


3.46544 


0.9429 


6.103 


[1.4+100] M© 


3.46497 


0.9429 


6.101 


ISCO 


3.46410 


0.9428 


6.000 


[10+500] M @ 


3.46499 


0.9429 


6.084 


[1.4+500] M© 


3.46500 


0.9429 


6.083 



TABLE II: Dimensionless values for the energy E and angular momentum i z as defined in Eq. Q at the point of transition 
?trans- The values for the energy and angular momentum at ISCO have been included for reference. 



By construction, the plunge waveform is a good description of the waveform all the way to the event horizon. 
However, we shall attach a set of quasinormal RD modes (QNMs) at the effective light-ring as in the transition model 
case. Following Buonanno, et. al. [3l|, we attach the RD modes at the time when the orbital frequency (155)) peaks. 
The frequency of these ringdown modes depends on the mass and spin of the newly formed Kerr black hole after 



merger. We use the following EOB-bascd fit for the final mass and spin [31 



M f /M 
a f /M f = q f 



1 + (^8/9 - 1)7? - 0.498(±0.027) t? 2 , 
v^t? - 2.900(±0.065) t? 2 . 



(74) 
(75) 



We note that this one-parameter fit for the final spin of the post-merger Kerr BH, Eq. (|75l) . differs from the zero-spin 
limit of the spin-dependent expression used for the "transition model" , Eq. ([4"Tj) . The two expressions render similar 
results, but we will use Eq. ([751 to estimate the final spin of the post-merger Kerr BH in our EOB waveform model. 
The reason for this is two-fold, firstly it ensures that we use a consistent EOB framework for this second model and 
secondly, it is the conservative approach. We will use the EOB model both as a consistency check of the transition 
model and to assess the level of uncertainty in our results. The latter is best accomplished by using this alternative 
fit to the final spin to make the waveform as different as possible within results available in the literature. We have 
verified that the waveform models are not particularly sensitive to this choice, and the level of consistency we quote 
later encodes this uncertainty. The matching onto QNMs is obtained as before by imposing continuity of the EOB 
waveform, Eqs. ([72p and ([73"[). and all the higher order time derivatives that are necessary to determine the amplitudes 
and phases of the leading RD tone, two overtones and their respective twin modes. The strategy to carry out this 
matching was described in Section [V] 
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VII. DYNAMICS OF THE WAVEFORM MODELS 



We shall now combine the ingredients of the preceding sections to discuss the orbital evolution and waveforms 
for our twelve test systems, which have mass ratios 77 = 0.0192 (1OM + 500M Q ), r] = 0.0028 (1.4M + 500M Q ), 
r/ = 0.0826 (10M Q + 100M Q ) and r] = 0.0136 (1.4M + 100M Q ), and the three different spin parameters, namely, 
q = 0, 0.3, 0.9. For the zero spin case, we will compare the EOB model to the transition model as a consistency check. 
The final ingredient we need to include in the model is the response function and the noise model for the ET detector. 



A. Implementation of the response function 

The ET response may be written as 

K{t) = ±[F+it)h + it)+FZit)h*it)], (76) 

where a = I, II refers to the two independent right-angle Michelson-like detectors. Any number of coplanar and 
colocated detectors with uncorrelated noise have an equivalent GW response to that of two right-angle interferometers 
offset by 45° to one another, and it is these putative detectors that we label as I and II. The functions h + ,x (t) are 
the two independent polarizations of the gravitational waveform. The antenna pattern functions F+ x are given by 



(77) 



(78) 

The various angles in the previous expressions represent, (a) the source's sky location in a detector based coordinate 
system (0,0), and (b) the polarization angle of the wavefront (?/>). These can be re-written in a fixed, ecliptic- 
based coordinate system. If we denote the source co-latitude and azimuth angles in the ecliptic coordinate system 
by i9s,(f>s), and the direction of the IMBH's spin a by then we can use the expressions (5), (7b) and 

(64a)-(64c) of Apostolatos et. al., [6(| to determine 9 it), tpit), ^>(i). 

If <I>(i) denotes the phase of the waveform at the centre of the Earth, there is also a shift in the waveform phase 
due to the difference in location of the detector. This can be included as a phase shift [6l[ 

$(*) -> $(*) 4 2^i?sin<9 s cos[27r(t/T) - <f> s ], (79) 

where R = i?©/c = 0.02125 s, and d<fi/dt is the azimuthal velocity of the orbit. This term is different for different 
detectors in a network, and encodes the time delay information that allows source triangulation using the network. 
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cos(20) cos(2-0). 



B. Noise model 



The signal-to-noise ratio for a given waveform is determined by an integral in Fourier space, weighted by the power 
spectral density (PSD) of the detector. For a monochromatic source, we can approximate this integral as a time 
domain integral by defining a noise-weighted waveform 

Kit) ^ -tML, /(,) = I§. (80) 

7r at 



s h (fit)) 

Where the PSD, Sh(f), is evaluated at the instantaneous gravitational wave frequency, fit). This is a good ap- 
proximation during the inspiral, during which the frequency is given by Eq. ([8]), and for the ringdown, during which 
the PSD factor takes the form (/(£)) = Shifimn), where fi mn — uj£ mn /2n [38]. Since the distorted Kerr BH 
will predominantly be emitting GWs at the leading mode I = m = 2,n = 0, we set Sh if it)) — Sh (/220) during 
the RD phase. In the merger and plunge phase, this approximation is not appropriate, but this phase is short and 
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so we use the same approximation, taking the instantaneous frequency from (|56|) in the EOB case and (|35|) in the 
transition model. Treated in this way, there is a small jump in the frequency used to evaluate the PSD when the 
orbit passes from plunge to ringdown, so we use interpolation to ensure a smooth transition. In the regime where 
this approximation is used, the PSD is quite flat (cf. Figure [T]) so we do not expect significant errors from using this 
approach. Furthermore, we have checked that this treatment of the detector noise generates results that are consistent 
with results computed directly in the frequency domain. 

C. Sample waveforms 

We shall now put at work the machinery developed in Sections IIIIHVII to generate a few sample noise-weighted 
waveforms. In Figure [2] we show the waveform for two sample systems. The waveforms in Figure [2] have the expected 
form. We see a gradual chirping signal with increasing frequency and amplitude which peaks at the merger and is 
damped exponentially afterwards. The ringdown radiation is weaker for the lower mass inspirals, m = 1.4M©, as we 
would expect since the energy released scales as r/ 2 . 



60 | 1 , , 1 1 1 , 1 8 




, 1 1 1 1 1 1 1 1 _i l 1 1 1 1 1 1 

7.5 7.55 7.6 7.65 7.7 7.75 7.8 7.85 7.9 54.6 54.65 54.7 54.75 54.8 54.85 54.9 

Time [seconds] Time [seconds] 



FIG. 2: Complete gravitational waveforms in the interferometer T for COs of masses 10 Mq (left panel) and IAMq (right 
panel), orbiting around a 5OOA/0 BH with spin parameter q — 0.3. The gravitational waveform shows the last stage of inspiral, 
plus the transition, plunge and RD phases. The various extrinsic parameters were chosen randomly. 




Time [seconds] 

FIG. 3: We show gravitational waveforms emitted during the final stage of inspiral, merger and ringdown for a 10 Mq CO 
orbiting around a 5OOA/0 BH, as computed using the EOB model and the transition waveform with q set to 0. As for Figure [5] 
the various extrinsic parameters were chosen randomly. 
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In Figure [3] we compare the EOB and transition models for the waveform. Since both models were designed to 
yield accurate results in the test-mass limit, we expect a good agreement for small 77, but perhaps a larger deviation 
for large r\. In fact, the models agree well even for the larger value of 77. The phase is in very good agreement, which 
is important for matched filtering, and the amplitude agrees to better than 10%. This provides confidence in our 
results, and the difference provides a guide to the overall level of uncertainty that we might expect in the SNRs which 
we quote in the next section. 

In the following section we present SNRs for these twelve different binary systems. The full parameter space of 
the waveform is ten dimensional and these parameters are defined in Table HTT1 For the SNR calculations, we fix the 
intrinsic parameters (the first four parameters of Table IIII[) and the distance to the source and then run a Monte 
Carlo over the values of the remaining five extrinsic parameters. 



lnm 


mass of CO 


InM 


mass of SMBH 


1 


magnitude of (specific) spin angular momentum of SMBH 


to 


time at which orbital frequency sweeps through a reference value 


00 


initial phase of CO orbit 


03 


source sky colatitude in an ecliptic-based system 


<t>S 


source sky azimuth in an ecliptic-based system 


6 K 


direction of SMBH spin (colatitude) 




direction of SMBH spin (azimuth) 


InD 


distance to source 



TABLE III: This table shows the physical meaning of the parameters used in our model. The various angles (9g,(f>g) and {Qki^k) are 
defined in an ecliptic-based coordinate system. 



VIII. SNRS FOR SPINNING AND NON-SPINNING BINARIES 



We start with a brief review of the framework of signal analysis. A GW detector can be thought of as a linear 
system whose input is a GW and whose output is a time series. This time series is a combination of both instrumental 
noise and a true GW signal. For ET, the output of each of the equivalent Michelson detectors can be represented as 

s a (t) = h a (t) +n a (t), a = 1,11. (81) 

The detection problem is to distinguish h a (t) from n a (t). We can define 



s(t)K(t)dt, (82) 

where K{t) is a filter function. Given the fact that we have a model for the signal h(t), we want to find the filter 
function that maximizes the SNR for such a GW signal. This optimal filter function is usually known as Wiener filter. 
There is a natural inner product on the vector space of signals, which for any two signals Paft), Qa(t), takes the form 

poo 

(p|q) = 2£/ {p* a (f)q a (f)+PM)?M)}/S h (f)df, (83) 

and in terms of this inner product, the maximal SNR obtained with the Weiner optimal filter is 

^W)l = ( s|h) (84) 

where the expected waveform h(t; 6) depends on parameters 6 = {9\, On}. 

If we consider white noise, so that Sh(f) is a constant, we can use Parseval's theorem to rewrite the inner product 
(|83| as 25 1 " 1 J2 a J-00 P<* (*)?<*(*) and the optimal value of the SNR becomes 



SNR 2 = 2 j 



*LSO 

h 2 a (t)dt. (85) 
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We can use this expression for the SNR if we use noise-weighted time-domain waveforms, as defined in Eq. (|80l) . 
We have used Eq. (|551 to compute SNRs for twelve sample binary systems at a fixed distance D — 6.63481 Gpc, 
which corresponds to the luminosity distance to redshift z = 1. As mentioned earlier, the triangular design of ET, 
consisting of three interferometers with 60 degree opening angles, generates a response equivalent to two co-located 
interferometers with 90 degree opening angles, but rotated 45 degrees with respect to each other, and with SNR that 
is a factor 3/(2\/2) ~ 1-06 higher. We have not included this factor in the quoted SNRs, given the uncertainties in 
the ET design that exist at this stage. The distribution of SNRs over random choices of the extrinsic parameters are 
summarised in Figures 0HZ1 while the statistics of the SNR distributions are given in Tables IIVI and [V] 




1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.8 1 1.2 1.4 1.6 1.8 2 

log(SNR) log(SNR) 



FIG. 4: The left panel shows the distributions of SNRs in a Monte Carlo simulation over extrinsic parameters for a 10 Mq + 
IOOMq binary with IMBH spin parameter q = 0.9, 0.3. The right panel shows the SNR distribution for the same binary 
system using the transition waveform model in the q — limit along with its EOB counterpart. The horizontal axis is the 
logarithm to base ten of the SNR. Note that the SNR distributions have been computed at a fixed distance D = 6.63481 Gpc, 
or, equivalently, at a fixed redshift z — 1. The top horizontal axis mirrors the values used in the bottom horizontal one, but in 
normal units. 




log(SNR) log(SNR) 

FIG. 5: As Figureg] but now for binaries with masses 1.4 Mq + 100M Q . The left panel shows SNR distributions for IMBH 
spin q = 0.3 and q = 0.9, while the right panel shows SNR distributions for IMBH spin q = computed using both the EOB 
and transition models. 

We can better visualize the results shown on Tables IIVI and [V] by plotting the mean of the SNR distributions as a 
function of the spin parameter q, which is shown in Figure [SJ We see that rapidly spinning binaries, q ~ 0.9, with 
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0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 0.4 0.6 0.8 1 1.2 1.4 1.6 



log(SNR) log(SNR) 

FIG. 6: The left panel shows the SNR distributions for a 10 Mq CO orbiting around a 500Afo BH with spin parameters 
q — 0.9, 0.3, respectively. The right panel shows the SNR distributions for the same binary system with q = using both 
waveform models. The panels show the distribution of the logarithm to base ten of the SNR. The SNR distributions have been 
computed at a fixed distance D = 6.63481 Gpc. The top horizontal axis mirrors the values used in the bottom horizontal one, 
but in normal units. 




log(SNR) log(SNR) 



FIG. 7: As Figured! but now for binaries with masses 1.4 Mq + 500M Q . The left panel shows SNR distributions for IMBH 
spin q = 0.3 and q = 0.9, while the right panel shows SNR distributions for IMBH spin q = computed using both the EOB 
and transition models. 



large mass ratios, rj ~ 0.08, will be relatively loud. At a fixed SNR detection threshold of 10, we conclude that such 
sources would be seen to distances D ;> 6.6 Gpc. In contrast, slowly rotating binaries with q ~ 0.3, and small mass 
ratios, r\ ~ 0.003, will only be visible at distances D <; 6.6 Gpc. 

In order to understand the implication of these results, we can use the data summarized on Tables IIVI and [V] to 
estimate the number of events per year that could be detected by the ET. To do so we will follow the procedure outlined 
in Section 3.3 of Q. The basic idea is the following: fix an SNR detection threshold, pthrcsh (we use pthresh = 10 in 
the following), then estimate the luminosity distance, Dl{z), at which a given source can be detected, i.e., have this 
particular SNR, as Dl = p(6.63481)/pthreshj where p(6. 63481) is the SNR of the source at a distance of 6.63481Gpc 
given in the Tables. We then use the concordance cosmology to convert this luminosity distance estimate into a 
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Stats 


m = 10M© 


m = 1.4M© 


q = 0.9 


q = 0.3 


q = 


EOB 


q = 0.9 


q = 0.3 


q = 


EOB 


Mean 


74.645 


46.774 


37.844 


38.282 


24.099 


14.388 


11.272 


11.967 


St. Dev. 


31.030 


20.665 


12.401 


14.450 


9.113 


5.149 


4.101 


4.103 


L. Qt. 


55.081 


33.266 


27.606 


27.669 


18.113 


10.740 


8.337 


8.933 


U. Qt. 


103.276 


65.163 


52.000 


51.642 


32.584 


19.454 


15.205 


16.144 


Med. 


73.451 


46.026 


37.844 


38.107 


24.266 


14.289 


11.041 


11.995 



TABLE IV: Summary statistics of the SNR distributions for binary systems with a central IMBH of mass M = IOOMq, and 
various choices for IMBH spin, q, and CO mass, m. We show the mean, standard deviation, median and quartiles of the 
distribution of the SNR for each system. The SNR distributions have been computed at a fixed distance D — 6.63481 Gpc. 
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m = 10M Q 


m = 1.4M© 


q = 0.9 


q = 0.3 


q = 


EOB 


q = 0.9 


q = 0.3 


5 = 


EOB 


Mean 


55.463 


18.408 


12.853 


13.677 


15.417 


3.908 


2.570 


2.729 


St. Dev. 


21.337 


6.723 


4.271 


4.716 


5.571 


1.256 


1.102 


1.113 


L.Qt. 


39.811 


13.583 


9.484 


9.908 


11.508 


2.884 


1.932 


2.084 


U.Qt. 


74.645 


25.645 


17.579 


19.099 


20.845 


5.272 


3.327 


3.707 


Med. 


55.590 


18.408 


12.764 


13.836 


15.668 


3.882 


2.636 


2.723 



TABLE V: As Table HVl but now for binary systems with a central IMBH of mass M = 5OOM . 




bj , , , i i i i i , , , i i i i i 
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FIG. 8: The plot shows the mean value of the SNR distribution as a function of the spin parameter q. From top to bottom, 
the various lines correspond to binaries of mass-ratio 77, namely, solid line 77 — 0.0826 ([10+100]Mq), dashed lined 77 — 0.0192 
([10+500]M Q ), dash-dot line 77 = 0.0136 ([1.4+1OO]M ) and dotted line 77 = 0.0028 ( [1. 4+500] M ). 



redshift estimate z by inverting the following expression, 

dz' 



D L (z) = D H (l + z)^ 



\n M {i + z'Y + n A } 1/2 



(86) 



We will assume a flat universe ilk — 0, and use J7m = 0.27, J7a = 0.73, Ho = 72 km s 1 Mpc , and Dh = c/Hq rj 
4170 Mpc. We can then compute the comoving volume, V c , within which the source can be detected using 62], 

v = ^p%[r dz' I 3 

3 \J in M (i + z>r + n A } 1/2 1 

In a previous study on the possible detection of IMRIs of COs into IMBHs with Advanced LIGO [63], it was found 
that binary tightening via three-body interactions was the dominant mechanism that led to the formation of IMRIs. 
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In this mechanism, the merger time for an IMRI can be estimated as the sum of the hardening timescale, Thard. 
and the gravitational-wave merger timescale, Tqwj he, Emerge = ^harden + ^GW 0; where 




in which a is the semimajor axis of the binary, n is the number density of stars in a globular cluster, a is the velocity 
dispersion, and m* stands for the mass of stars that interact with the binary. In practice, we set n, a and m* to their 
fiducial values, and minimize T merge = Th ar den + 7gw over a to estimate the CO-IMBH coalescence rate. The rate at 
which IMRIs occur per globular cluster can then be approximated by l/T mergo . 

The SNRs that we calculate are the SNRs for systems with redshifted masses, M z = M(l + z), m z = m(l + z), 
equal to those that we have specified. The maximum detectable redshift estimated from the SNR then tells us the 
intrinsic source-frame masses of the system whose range we have computed. In the case of ET, which can detect 
sources out to cosmological redshifts, z ;> 1, these intrinsic masses may not correspond to astrophysically interesting 
systems. However, by considering a range of redshifted masses and computing for each one an event rate under the 
assumption that all IMRI systems were of that particular intrinsic type, we can still obtain a rough estimate of the 
event rate. 

Following Q, we assume that 10% of clusters form an IMBH and are sufficiently dense to be hosts to an IMRI and 
we assume that such globular clusters have a fixed comoving density of ~ 0.3Mpc -3 . The rate of detectable events for 
a particular type of system can then be estimated as ~ 0.3(K/Mpc 3 )/[T mergc (l + z)} [7]. Table IVT1 presents estimates, 
for each of our fiducial systems, of the maximum detectable redshift, the intrinsic masses that the source represents 
at that redshift, the corresponding merger time through binary hardening, the comoving volume within that redshift 
and the IMRI event rate assuming all IMRI sources were of that type. Note that the latter assumption means that the 
entries in this table are not independent of each other, i.e., the total number of events is not given by the sum of the 
last column, but is somewhere in the range of rates tabulated. Even though the astrophysical properties of IMBHs, 
e.g., their mass and spin distributions, are currently very uncertain, we can still draw conservative predictions from 
Table IVII These estimates compare fairly well with those of [tJ , but are a bit bigger since we use a more accurate 
waveform model. We see that IMRIs could be seen at redshifts z ~ 1-6, depending on the mass of the IMBHs that 
exist, and there could be as many as a few hundred systems observed. If IMRIs tended to be 2M + 100M© systems, 
they could be seen out to z ~ 4 and we'd expect a few hundred events, but if IMRIs tended to be IMq + 400Af Q 
systems, we would only see ~ 10 events out to z ~ 0.3. The greatest uncertainty in these figures comes from the 
unknown number of IMBHs that exist in the Universe, and these uncertainties are not folded into the numbers in 
Table IVll If IMBHs are rare, then the IMRI rate could be orders of magnitude lower. 
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1.1 
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TABLE VI: "3 ET detector network" average range, corresponding redshift, source-frame masses, merger timescale, comoving 
volume within range, and detectable event rate for several combinations of plausible redshifted CO and IMBH masses. 



As discussed earlier, there is currently some uncertainty about the low-frequency sensitivity that ET will achieve. 
The preceding results assumed a frequency cut-off at 5Hz, but it is informative to consider how the SNR changes if 
this frequency cut-off is pushed down to 3Hz or even 1Hz. In Table [VlT] we show how the SNR changes as a function 
of the cut-off frequency, for each of the binary systems. In each case we have chosen a particular random, but fixed, 
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set of extrinsic parameters that is representative of the events that lie around the peak of the SNR distribution for 
the binary. 



Binary 


q = 0.9 


g = 0.3 


q = 


f = 5Hz 


f = 3Hz 


f = 1Hz 


f = 5Hz 


f = 3Hz 


f = 1Hz 


f = 5Hz 


f = 3Hz 


f = 1Hz 


[10 + 100] M© 


72.111 


79.799 


86.896 


52.240 


56.885 


63.680 


42.855 


45.186 


48.306 


[10 + 500] M© 


55.286 


58.325 


63.293 


19.364 


20.701 


22.961 


15.276 


15.959 


16.331 


[1.4 + 100] Mg 


20.941 


21.878 


23.714 


16.711 


17.418 


18.030 


11.351 


11.952 


12.123 


[1.4 + 500] M 


13.274 


13.804 


13.932 


4.207 


4.256 


4.613 


2.610 


2.792 


3.029 



TABLE VII: Summary of how the SNR of a source changes as the low-frequency sensitivity cut-off changes from 5Hz to 3Hz to 
1Hz. For each binary, the extrinsic parameters have been chosen randomly, but kept fixed as the cut-off frequency was varied. 

We see that if ET achieves frequency sensitivity down to 1Hz, we will have somewhat greater sensitivity to systems 
containing more massive compact objects. This will facilitate the extraction of these signals from the data, and open 
up the possibility of detecting these sources at higher redshift. We have checked the reliability of the figures quoted 
in Table IVTlI by running a full Monte Carlo for the particular system [10 + 100] Mq, with spin parameter q = 0.3, and 
frequency cut-off at 1Hz. Figure |9] presents these results, and confirms that Table IVlTl is a fair representation of what 
may be achieved by ET at lower frequencies. 



15.8 | 39,8 | 100 




log(SNR) 

FIG. 9: The panel shows the SNR distribution for a 10 M© + 100M© binary with IMBH spin parameter q — 0.3, at a frequency 
cut-off of 1Hz. The horizontal axis is the logarithm to base ten of the SNR. Note that the SNR distribution has been computed 
at a fixed distance D — 6.63481 Gpc, or, equivalently, at a fixed redshift z = 1. The top horizontal axis mirrors the values used 
in the bottom horizontal one, but in normal units. The statistics of the system are summarized at the top-left of the panel. 

We can use the figures quoted in Table IVIII to explore the redshift at which the loudest sources could be seen if 
ET achieved a frequency cut-off of 1Hz. From Table l*VTl we already know that the systems with redshifted masses, 
m z + M z , of [10 + 1OO]M , [10 + 500]A/©, with q = 0.9, could be seen up to redshift z ~ 5, z ~ 4, respectively. 
Using Eq. ([86]). we see that, with a cut-off at 1Hz, these same systems could be detected up to redshift z ~ 6, z ~ 5, 
respectively. Table IVlTTl summarizes these results for the loudest sources from Table fVTl The event rate estimate only 
changes appreciably for the q = 0.3 case (since the intrinsic masses for the systems are also changing), but as q = 0.3 
may be a good estimate for the typical spin of an IMRI, there could be a significant scientific gain from pushing the 
lower frequency cut-off to 1Hz. A more systematic study, which fixes the intrinsic masses of the events as opposed to 
the redshifted masses, is needed to fully explore the implications of the cut-off on the expected IMRI detection rate. 

There is another important consequence of a lower low-frequency cut-off. As discussed previously, the binary 
systems we have considered in this paper are very short lived. However, by pushing the seismic wall down to 1Hz 
these binaries now stay in the sensitivity band of the detector for longer. Therefore, we should be able to obtain 
better parameter estimation accuracies, in particular for the extrinsic parameters. Table HXl gives the time that each 
of these binaries is in band as a function of the frequency cut-off. 

Table llXl shows that by pushing the frequency cut-off to 1Hz will boost the time spent in band of the shortest-lived 
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M z /M e 


m z /M 


q 


D/Gpc 


z 


M/Mq 


m/M 


Tmcrgc / yr 


V c /Mpc 3 


Events/yr 


100 


10 


0.9 


57.72 


5.9 


14.5 


1.4 


5.81 x 10 s 


2.47 x 10 12 


187 


100 


10 


0.3 


42.46 


4.6 


17.8 


1.8 


4.47 x 10 s 


1.92 x 10 12 


202 


500 


10 


0.9 


41.79 


4.5 


90.9 


1.8 


2.62 x 10 s 


1.87 x 10 12 


385 


500 


10 


0.3 


15.26 


2.2 


156.2 


3.1 


1.70 x 10 s 


7.08 x 10 11 


346 



TABLE VIII: "3 ET detector network" average range, corresponding redshift, source-frame masses, merger timescale, comoving 
volume within range, and detectable event rate for several combinations of plausible redshifted CO and IMBH masses. The 
cutt-off frequency has been set at 1Hz. 



Binary 


q = 0.9 


<7 = 0.3 


q = 


f = 5Hz 


f = 3Hz 


f = 1Hz 


f = 5Hz 


f = 3Hz 


f = 1Hz 


f = 5Hz 


f = 3Hz 


f = 1Hz 


[10 + 100] M 


45.0 


169.2 


3093.7 


42.1 


161.5 


3048.5 


38.9 


157.3 


3025.8 


[10 + 500] M 


16.3 


61.2 


1099.3 


8.3 


43.9 


1010.1 


5.4 


35.6 


963.2 


[1.4 + 100] M e 


319.9 


1209.0 


22089.1 


291.0 


1152.4 


21776.9 


275.9 


1123.7 


21614.1 


[1.4 + 500] M Q 


112.5 


436.1 


7851.8 


55.6 


311.2 


7221.9 


34.0 


252.5 


6879.9 



TABLE IX: Summary of how the time spent in band of a source changes as the low- frequency sensitivity cut-off changes from 
5Hz to 3Hz to 1Hz. The time is reported in seconds. For each binary, the extrinsic parameters have been chosen randomly, 
but kept fixed as the cut-off frequency was varied. 

events by a factor of ~ 200. Binaries of [1.4 + 100] Af© could spend up to 6 hours in band. In the second paper of 
this series we will show that even with a 5Hz cut-off we can pinpoint the location of these particular sources in the 
sky and constrain their luminosity distances with a precision of ~ 10%, at SNR of 30, using a 3 ET detector network. 
Pushing the seismic limit down to 1Hz will further improve the parameter estimation accuracies. 

All the results obtained up to this point have assumed that we will be able to detect IMRIs using a 3 ET detector 
network. Hence, our previous results may be considered as upper limits for the various quantities we have quoted 
— event rates per year, SNRs, horizon distances etc. We shall now relax this assumption and explore more modest 
scenarios, the configurations C1-C5 outlined earlier. To recap, these are, CI: one ET at the geographic location of 
Virgo; C2: as configuration CI plus a right-angle detector at the location of LIGO Livingston; C3: as configuration 
CI plus another ET at the location of LIGO Livingston; C4: as configuration C2 plus another right-angle detector 
in Perth; and C5: as configuration C3 plus another ET in Perth. Note that configuration C5 corresponds to the 3ET 
detector network used in our analyses so far. We will quote results for IMBHs with spin parameter q ~ 0.3 only, since 
as argued earlier this could be a reasonable fiducial value for IMBH spin. We quote SNRs for IMRIs into a central 
black hole of mass M — IOOMq and spin q = 0.3 for five different configurations in Table Ixl and corresponding results 
for IMRIs with M = 500M Q and q = 0.3 in Table EH 



Stats 


m = 10M Q 


m = 1.4M Q 


CI 


C2 


C3 


C4 


C5 


CI 


C2 


C3 


C4 


C5 


Mean 


26.242 


34.119 


40.179 


40.458 


46.774 


8.222 


10.740 


12.589 


12.823 


14.388 


St. Dev. 


14.785 


16.623 


18.912 


19.015 


20.665 


3.113 


3.485 


4.102 


4.408 


5.149 


L.Qt. 


18.239 


24.660 


28.708 


29.242 


33.266 


5.754 


8.054 


9.333 


9.311 


10.740 


U. Qt. 


38.282 


48.195 


57.148 


55.719 


65.163 


11.169 


14.928 


16.711 


16.634 


19.454 


Med. 


26.424 


33.884 


39.719 


40.365 


46.026 


8.017 


10.666 


12.162 


12.445 


14.289 



TABLE X: As Table IIVI but for binary systems with a central IMBH of mass M — 1OOM and spin parameter q = 0.3 and 
assuming four additional configurations for the detector network, C1-C5 as described in Section Til Al Configuration C5 is the 
network of three ETs which has been used for all results elsewhere in this paper. 

From Tabl es IXl and IXI I we learn that the SNR corresponding to a 3 ET detector network is, roughly speaking, a factor 
of y/3 and a/3/2 greater than the SNR associated with a single ET and a 2 ET network, respectively. We may expect 
this scaling since, despite the different locations of the detectors, we would expect the SNR to scale approximately as 
the square root of the number of detectors. We also notice that the SNR associated with configuration C4 is similar 
to the SNR of the 3 ET network. This configuration is somewhat less ambitious than a 3ET network and hence might 
be more likely to be realised in the future. We will show in the second article of this series that for slowly rotating 
IMBHs, such a network will also suffice for determining the luminosity distance to an event with an accuracy of ~ 12% 
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m = 10M© 


m = 1.4M 


CI 


C2 


C3 


C4 
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CI 


C2 


C3 


C4 


C5 


Mean 


9.528 


12.677 


14.825 


14.894 


18.408 


2.084 


2.685 


3.141 


3.177 


3.908 


St.Dev. 


4.432 


4.866 


5.551 


5.346 


6.723 


1.102 


1.117 


1.164 


1.129 


1.256 


L.Qt. 


6.683 


8.974 


10.328 


12.677 


13.583 


1.517 


2.032 


2.333 


2.382 


2.884 


U.Qt. 


13.868 


17.742 


20.370 


20.045 


25.645 


2.844 


3.589 


4.188 


4.207 


5.272 


Med. 


9.638 


12.853 


14.521 


14.689 


18.408 


2.109 


2.704 


3.090 


3.126 


3.882 



TABLE XI: As Table IXl but for binary systems with a central IMBH of mass M — 5OOM0 and spin parameter q = 0.3. 

at a source SNR of 30. If the luminosity distance is converted into a redshift using the concordance cosmology at that 
time, then the resulting redshift error will be comparable to the distance error. Furthermore, configuration C4 will 
also allow us to determine the masses of the components of the system to high accuracy. 

The above discussion indicates that ET, through the detection of gravitational waves emitted by IMRIs of COs 
into IMBHs, may shed some light on the astrophysical properties of IMBHs, e.g., their mass and spin distributions, 
and to find out whether dynamical interactions in dense stellar systems do indeed lead to the formation of IMBHs. 
Additionally, the mergers detected by ET will be complementary to mergers between heavier BHs that will be seen 
by space-based detectors such as LISA, ALIA or DECIGO. These various detectors will provide a measurement of 
the rate of black hole mergers in various mass ranges, which will be useful to place constraints on models of black 
hole growth 0- 

IX. CONCLUSIONS 

We have developed waveform models for IMRIs of COs on circular, equatorial orbits into central IMBHs of arbitrary 
spin. One waveform model, which is valid for inspirals into IMBHs of arbitrary spin, uses the transition-to-plunge 
scheme of Ori and Thorne [l| to smoothly match the inspiral onto a plunge waveform and ringdown. The second 
approach, at present valid only for non-spinning IMBHs, uses the effective-one-body formalism to match a merger and 
ringdown onto the inspiral waveform. We have shown that the two distinct waveform families are in good agreement 
in the q = limit, particularly in the small rj regime. The agreement in phasing all the way from inspiral to ringdown 
is particularly good and this is important as a good phase model will be crucial for the detection of these systems via 
matched filtering and to extract parameter information from the detector measurements. 

We have used these waveform models to compute estimates for the SNR that would be obtained in ET for various 
binary systems. We computed SNRs for a 3 ET detector network for twelve different binaries using the transition 
model and used the EOB model to cross-check the results for the non-spinning IMBHs. We found that the two models 
made predictions that were consistent to about ten percent. Assuming that ET has a low- frequency sensitivity cut- 
off at 5Hz, we found that at a redshift of z = 1, typical SNRs for IMRI systems with masses 1.4Mq+100Mq, 
1OM q +1OOM , 1.4M q +500M q and 10M Q +500M Q will be in the range ~ 10-25, ~ 40-80, ~ 3-15 and - 10-60 
respectively. Using the SNR distributions as input data, we estimated the horizon distance at which these various 
sources could be seen. This suggested that ET could detect as many as several hundred of these systems, up to a 
redshift z <J 5, although the exact number will depend on the intrinsic distribution of masses and spins for the IMRI 
systems. If the ET sensitivity extends down to 1Hz, we found that the same systems could be detected up to redshift 

We have also explored more modest network configurations consisting of 1 ET only, 1 ET plus 1 right-angle 
detector, 2 ETs and 1 ET plus 2 right-angle detectors. Using these configurations we computed SNR distributions 
for the same four combinations of source masses 1.4M +1OOM , 10M Q +100M Q , 1.4M +5OOM Q , 10M Q +500M Q , 
but fixed IMBH spin parameter q = 0.3. We chose this spin since it is a reasonable estimate of the spin parameter 
of an IMBH that has grown primarily through a series of minor mergers. A network consisting of one ET and two 
right-angle interferometers will have almost as great a sensitivity to IMRIs as the highly-ambitious 3-ET network, 
and this result will be important when a third-generation detector network is planned. 

We shall see in the second paper of this series that a configuration consisting of 1 ET only will not be enough to 
effectively constrain the distance to the source, let alone to pinpoint the location of the source in the sky. However, 
the configuration consisting of 1ET and 2 right-angle detectors could determine the luminosity distance of a source 
to an accuracy of ~ 15% at SNR of 30. Since ET will also be able to determine the masses of a binary to great 
accuracy, we should be able to estimate with confidence the masses and the redshift at which two COs merge. This 
information will be important for understanding the astrophysical properties and history of IMBHs. IMRI detections 
will probe the existence and properties of IMBHs, and their number density over cosmic history. This in turn will 
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tell us how these objects form and evolve. Since IMRIs will primarily be observed from globular clusters, IMRI 
observations will tell us about the efficiency of formation of IMBHs in cluster environments and the number density of 
cluster IMBHs 0] • This in turn will allow us to estimate the rate at which IMBH-IMBH binaries could form through 
the globular cluster channel. ET might also detect IMBH-IMBH binaries that arise from primordial IMBHs |64| . 
Understanding the properties of globular cluster IMBHs through the IMRI channel could therefore help to identify 
candidate primordial IMBHs, which will have important consequences for hierarchical models of structure growth. 
IMRIs could also be used to test gravitational physics in the strong field and map the structure of spacetime outside 
IMBHs [7] . The results in this paper indicate that ET might detect as many as several hundred IMRI events which 
could be used to extract this physics. These results also lay a solid foundation for paper II in this series, which will 
study the precision with which the parameters of these IMRI systems might be estimated using the ET detector. 
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